home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dlas2.z / dlas2
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAASSSS2222((((3333FFFF))))                                                            DDDDLLLLAAAASSSS2222((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLAS2 - compute the singular values of the 2-by-2 matrix  [ F G ]  [ 0 H
  10.      ]
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DLAS2( F, G, H, SSMIN, SSMAX )
  14.  
  15.          DOUBLE        PRECISION F, G, H, SSMAX, SSMIN
  16.  
  17. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  18.      DLAS2  computes the singular values of the 2-by-2 matrix
  19.         [  F   G  ]
  20.         [  0   H  ].  On return, SSMIN is the smaller singular value and SSMAX
  21.      is the larger singular value.
  22.  
  23.  
  24. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  25.      F       (input) DOUBLE PRECISION
  26.              The (1,1) element of the 2-by-2 matrix.
  27.  
  28.      G       (input) DOUBLE PRECISION
  29.              The (1,2) element of the 2-by-2 matrix.
  30.  
  31.      H       (input) DOUBLE PRECISION
  32.              The (2,2) element of the 2-by-2 matrix.
  33.  
  34.      SSMIN   (output) DOUBLE PRECISION
  35.              The smaller singular value.
  36.  
  37.      SSMAX   (output) DOUBLE PRECISION
  38.              The larger singular value.
  39.  
  40. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  41.      Barring over/underflow, all output quantities are correct to within a few
  42.      units in the last place (ulps), even in the absence of a guard digit in
  43.      addition/subtraction.
  44.  
  45.      In IEEE arithmetic, the code works correctly if one matrix element is
  46.      infinite.
  47.  
  48.      Overflow will not occur unless the largest singular value itself
  49.      overflows, or is within a few ulps of overflow. (On machines with partial
  50.      overflow, like the Cray, overflow may occur if the largest singular value
  51.      is within a factor of 2 of overflow.)
  52.  
  53.      Underflow is harmless if underflow is gradual. Otherwise, results may
  54.      correspond to a matrix modified by perturbations of size near the
  55.      underflow threshold.
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.